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Abstract 



We calculate the heavy quark potential from the magnetic current due to monopoles in 
four dimensional U{1) lattice gauge theory. The magnetic current is found from link 
angle configurations using the DeGrand-Toussaint identification method. The hnk angle 
configurations are generated in a cosine action simulation on a 24^ lattice. The magnetic 
current is resolved into large loops which wrap around the lattice and simple loops which 
do not. Wrapping loops are found only in the confined phase. It is shown that the long 
range part of the heavy quark potential, in particular the string tension, can be calculated 
solely from the large, wrapping loops of magnetic current. 



In this paper, we report new results on confinement via monopoles for U{1) lattice 
gauge theory in four dimensions. Our main result is that the confining part of the heavy 
quark potential, in particular the string tension, is determined solely by large loops of 
magnetic current. It has been established for some time that large loops which extend 
over the entire lattice are present only in the confined phase of the theory |[. Their 
presence can now be quantitatively tied to the string tension. Our work is carried out on 
a 24^ lattice, near the deconfining transition. 

The role of monopoles in U{1) lattice gauge theory is seen most clearly using the 
Villain [0] or periodic gaussian form of the U{1) theory. Under a dual transformation, 
the usual link angle description goes over into one involving an integer-valued magnetic 
current m^(x), defined on the links of the dual lattice [Q. The link angle path integral 
becomes a sum over all possible configurations of magnetic current. In this monopole 
representation, the system can be visualized as a plasma of magnetic monopoles moving 
on Euclidean world lines, interacting via photon exchange. 

In either representation, a Wilson loop calculation is needed to determine the heavy 
quark potential. In the link angle representation, a Wilson loop is specified by the expo- 
nential of a line integral : 

W{R,T) = (exp (^tJ2e,{x)J,{x)^^ , (1) 

where the integer-valued electric current is non-vanishing on the rectangular R x T 
loop contour, and < ■ >0 denotes the expected value taken over configurations of link 
angles 6^{x). 

In the monopole representation, the determination of a Wilson loop involves the ex- 
ponential of an area integral over a surface with the loop contour as its boundary 0, ^. 
The electric current is first expressed as the curl of a Dirac sheet variable 0; where 
dy denotes a discrete derivative. The sheet variable D^^ is nonunique. For \J^\ = 1, a 
specific choice is made by setting D^^ = 1 on the plaquettes of an (arbitrary) open surface 
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with boundary J^, and D^^ = on all other plaquettes. The area integral represents the 
dual flux set up by the magnetic current through this surface. To compute it, we define 
the magnetic vector potential 

y 

where v satisfies —d-dv{x — y) = 5x,y The field strength is given by F^jy = S^A^ — S^A^. 
In terms of D^^ and the dual field strength F*j^{x) = ^e^uafB^apix), the monopole repre- 
sentation of a Wilson loop is finally given by: 

W{R, T) = W,^,{R, T) ■ ^exp D,.{x)F;^{x)^ ^ , (3) 

where < • >m denotes the sum over configurations of magnetic current. The factor 
27r which appears in the exponent of Eq. arises from the Dirac condition on the 
product of electric and magnetic charge, and guarantees that the value of a Wilson loop 
is independent of the surface chosen to define D^^. The prefactor in Eq. (Q) describes one 
photon exchange between the quark and anti-quark: 



Wph{R, T) = exp ^-y J2 Jf^{x)'v{x - y)Ji,{y)^ , 



(4) 



where the electric coupling is related to the coupling f3v in the Villain action by 
= 1/ Pv- The factor Wph{R,T) contributes a purely perturbative Coulomb term to the 
heavy quark potential. 

Fortunately, the numerical evaluation of Wilson loops via Eq. does not require a 
direct simulation in terms of the magnetic current m^(a;). This is impractical owing to 
the long-range interactions generated by photon exchange between the monopole currents. 
DeGrand and Toussaint [0 showed how to locate monopoles directly in configurations of 
link angles. In their procedure, the plaquette angle 0^^{x) = d^Oy — 0^6^ is resolved into 
a fluctuating part 6fj_!^{x), and an integer- valued Dirac sheet variable m*^,(x): 

O^^i^ix) = 9^u{x) + 27rm*^(x), (5) 
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where rrf^^^x) = j^nuais^aisix) and O^^y^x) G (— 7r,7r). The magnetic current is then given 
by m^{x) = d^m^yi^x). This procedure allows only values of G [±2, ±1,0], whereas 
in principle all integer values are allowed. However, at values of the coupling near the 
deconfining transition, the values = ±1 are overwhelmingly dominant; even = ±2 
occurs only a small fraction of the time. Thus negligible error is caused by omitting higher 
values of m^. 

The monopole current can be used to calculate physical quantities, as well as merely 
to count monopoles. In our previous work in four dimensions [^, we used the current 
m^(x) to evaluate Wilson loops from Eq. (^). Similarly in three dimensions 0, we used 
the monopole density m(x) to evaluate Wilson loops using the c? = 3 analog of Eq. (|). In 
both cases the resulting heavy quark potential agreed with that extracted from the link 
angles and Eq. (|^), to within statistical errors. In the present work on a 24*^ lattice, we 
have again checked that potentials deduced directly from link angles and Eq. (|l|) agree 
with those obtained from the magnetic current and Eq. @. These calculations show that 
quantitative results on confinement can be obtained using topological objects. 

The derivation of Eq. (|]) as an exact formula is only possible for Villain's form of 
the U{1) theory. On the other hand, Wilson's cosine form of the U{1) action can 
be simulated much more efficiently. In our previous work we have shown that 

Villain action results can be extracted from a cosine action simulation, if a simple coupling 
constant mapping is used. More precisely, a simulation using the cosine action at coupling 
[3w is equivalent to a Villain action simulation at coupling with I3v related to j3w 



where Jq and Ji are modified Bessel functions. Eq. (||) determines the value of l//?y, and 
hence e^, which result from a cosine action simulation at a given value of (3w- The factor 



cosine action link angle configurations and the result used to calculate the second factor 
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by 10 




(6) 



Wph in Eq. (|^) is then completely determined. The magnetic current is identified from the 



of Eq. (^. The Wilson loops calculated in this manner using the cosine action and Eq. (^) 
differ from pure Villain action results by a harmless perimeter term. The R-dependent 
terms in the potentials agree within statistical errors P, |lT]. 

To summarize, our simulations use the cosine form of the U{1) action, identify m^{x) 
using the DeGrand-Toussaint method, and evaluate Wilson loops in the monopole repre- 



sentation. The link angle configurations were generated using a heatbath algorithm [12 



The calculation of A^{x) from mn{x) in Eq. was done using a four dimensional vec- 
torized FFT |[13| , [14| . The R x T rectangle lying in the Wilson loop plane was used 
as the defining surface for Df^,y. Magnetic current configurations were saved every 10 
sweeps. After Wilson loops were obtained from these configurations using Eq. (^, po- 
tentials were extracted using standard methods. The heavy quark potential V{K) was 
obtained from a straight line fit of lnW{R,T) vs. T, over an interval Tmin{R) to Tmax-, 
where Tmin{R) = R + 2 ioi R = 2,3, and R + 1 otherwise, while Tmax = 16. To determine 
the string tension a and Coulomb coupling a, the potentials were then fitted to a linear 
plus Coulomb form, V{R) = aR — a/ R+Vq, over the interval R = 2 to R = 7. The large 
string tension present in U{1) makes it difficult to work at larger values of R. Errors in 
physical quantities were estimated using both the jacknife method and binning the data 
into bins of various size. 

It is well established for U{1) that appreciable correlation lengths occur only in the 
immediate vicinity of the deconfining phase transition. The location of the phase tran- 
sition moves to larger values of jSw as the lattice size increases, in a manner roughly 
consistent with finite size scaling theory 0, Since only lattices of size up to 16^ 

were available in the published literature when we began our work, it was first necessary 
for us to locate the transition for a 24^ lattice. To do this, we performed a series of 
runs with various initial configurations for 1.0100 < f3w < 1-0120, and monitored the 
value of the 1x1 Wilson loop, W{1,1). For Pw > 1.0114, the system always reached 
a state with W{1, 1) ~ 0.65. For P^y < 1.0112, the system always reached a state with 
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W{1, 1) ~ 0.63. Subsequent analysis of the heavy quark potential showed these two states 
to be deconfined and confined, respectively. While we have not precisely located the de- 
confining phase transition, consistency with our results requires that the transition be in 
the interval 1.0112 < (3w < 1.0114 for a 24^ lattice. 

To avoid problems associated with long autocorrelation times that occur near the 
transition, we chose to use a run of 20,000 sweeps at Pw = 1.0103 for the results to be 
presented below. At this value of Pw, the correlation length is large enough to observe 
the beginnings of continuum behavior, but small enough to avoid problems with long 
autocorrelation times. The autocorrelation time r measured from the 1x1 Wilson loop 
was approximately 100 sweeps. In Fig.(l), we show the heavy quark potential determined 
from Eq. for (3]y = 1.0103 using 936 configurations of magnetic current. A linear- 
plus-Coulomb fit gave a string tension of a = 0.56(1), and a Coulomb coupling of a = 
0.30(2). The total number of links carrying magnetic current at this Pw was 98, 400(800). 
For comparison, we also show in Fig.(l), the potential determined from Eq. (|^) for the 
deconfined Pw value, Pw = 1.0114, where 400 configurations of magnetic current were 
analyzed. At this value of Pw, the string tension was statistically zero, while the Coulomb 
coupling was a = 0.22(4). The total number of links carrying magnetic current was 
40,000(300). 

We now turn to the resolution of the magnetic current into loops. For every other 
configuration of magnetic current, or every 20 lattice sweeps, magnetic current loops were 
individually identified and catalogued. The loop-finding algorithm proceeded by choosing 
a non-zero current link m^(xo) and following the current it carried through the lattice until 
a loop was completed by a return to the site xq. This process was carried out repeatedly 
from different starting points and ended when the entire configuration of current had 
been resolved into loops, with each current- carrying link belonging to a specific loop. The 
algorithm was deterministic: when looking for an outgoing current link at a particular 
lattice site, the direction ^ = I was chosen first, followed by yU = 2,3,4 . Intersections 
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of loops did occur (i. e. more than one outgoing link associated with a site), so the set 
of loops identified was not unique. However, self-intersections of loops were relatively 
rare, occurring with approximately the same probability as self-intersections of a purely 
random walk in = 4 [|16[ . 

Each loop identified as described above forms a closed path of links and automatically 
satisfies current conservation. Two different topologies are possible, depending on whether 
tracing out a loop also involves winding or wrapping around the lattice. A lattice with 
periodic boundary conditions is a hypertoroid, so topologically nontrivial loops which wrap 
around the lattice are permitted. To distinguish the two possibilities, the net current was 
measured for each loop: 

xGloop 

It is easy to show that only loops which wrap around the lattice have a non- vanishing A^, 
and that the components of A^ must be integer multiples of the lattice size along an axis; 

= "^^1 ■ N, for a cubic lattice of size N^. The integer is the net current of the loop 
in the /ith direction. For example, for a loop with A4 7^ 0, the sum of the charge density 
m4 over each spacial cube or "time slice" will equal the net charge n^. Likewise for other 
directions. While an individual loop can have a non-vanishing A^, a net current cannot 
actually occur on a finite lattice, so the sum of A^ over all loops vanishes identically. In 
our data, it was typical for a loop with non- vanishing A^ to be wrapped around the lattice 
several times in more than one direction. Values of \n^\ up to 10 were observed. In what 
follows we will use the terms "wrapped" for loops with non-vanishing A^ and "simple" 
for loops with vanishing A^. 

An indication that wrapped loops are crucial for confinement is that they are present 
only in the confined state and never are observed to occur in the deconfined state. As 
mentioned earlier, there are almost twice as many links carrying magnetic current in the 
confined phase. When the current is broken up into loops, it is found that this excess in 
the confined phase consists of a small number of wrapped loops. At (3 = 1.0103, out of the 
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total of 98,400(800) links carrying magnetic current, 51,000(400) are in wrapped loops, 
the remainder in simple loops. The average number of wrapped loops is only 4.6(1), so 
that the typical wrapped loop contains thousands of links. In contrast the average total 
number of simple loops is 6, 210(6), of which 3, 507(4) are in the form of elementary one- 
plaquette current loops composed of four links. The number of simple loops with a given 
number of links decreases rapidly as the number of links increases. Over 90% of the links 
in simple loops are included in loops with 60 links or less. 

Since wrapped loops occur only in the confined phase, it is natural to ask if they can 
explain the long range, confining part of the heavy quark potential. To investigate this, 
we computed the heavy quark potential again using Eq. (^, but for each configuration, 
including only the magnetic current from wrapped loops. The results are shown (omitting 
the photon factor Wph) in Fig. (2). A linear plus Coulomb fit to the resulting potential 
gave a string tension = 0.56(2), and a Coulomb term = 0.09(1). The string tension 
is within statistical errors of the value 0.58(2) found earlier from the heavy quark potential 
calculated using the full magnetic current. Next, we carried out a similar calculation using 
only the magnetic current from the simple loops. This produced the rather flat potential 
also shown in Fig. (2). A linear plus Coulomb fit to this potential gives zero string tension 
within statistical errors (cr^ = .0000(6)), and a Coulomb term = 0.06(1). The result 
of these two fits gives strong evidence that in the long distance region, there is a clean 
separation between the contributions of the two classes of loops. Only the large, wrapped 
loops containing thousands of links contribute to the confining part of the potential. This 
has been demonstrated here only within certain error bars, but it may well be an exact 
statement. 

In the fits described above for (3^ = 1.0103, the wrapped loops required a Coulomb 
term with coupling = 0.09(1), while the simple loops required a Coulomb term with 
coupling as = 0.06(1). In addition, there is a Coulomb term coming from the Wph factor 
in Eq. (^. Using Eq. (^ to evaluate /5y(1.0103), gives aph = 0.13 as the Coulomb coupling 
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arising from Wph- Simply adding the various terms , we obtain aph + as + a^; — 0.28(2), 
consistent with our previous result of 0.30(2) obtained with the full magnetic current. 

The results on the string tension and Coulomb coupling are consistent with additivity 
of the potential over the various contributions. In Fig. (3), we compare the potential 
determined from the full magnetic current and Wph (shown previously in Fig.(l)), with 
the potential obtained by summing the contributions from Wph, wrapped loops, and 
simple loops, plus a constant. A glance at Fig. (3) shows that the agreement is quite 
good. Additivity of the potential over the various contributions would imply that the 
contributions from wrapped and simple loops factor in the average over configurations. 
To check factorization, we performed a fit to the "potential" extracted from the ratio of 
Wilson loops assuming factorization to Wilson loops calculated with the full magnetic 
current: 

{WUR, T))^ ■ {Ws{R, T))^ I {W^{R, T) ■ Ws{R, T))^ . (7) 

The string tension and Coulomb coupling resulting from this were both zero to within 
statistical errors. This shows that factorization and therefore additivity of the potential is 
consistent with our data. This is not surprising at large R, where the R-dependence comes 
predominantly from the wrapped loops. However, in the small R region, both wrapped 
and simple loops produce Coulomb terms, and additivity is not expected to hold as an 
exact statement. Nevertheless it appears to be a good approximation and holds within 
the accuracy of our data. 

We have shown that confinement in t/(l) comes about through a particular component 
of the magnetic current, the large wrapping loops. The small, simple loops contribute 
only to the Coulombic part of the potential. While the latter is intuitively reasonable, 
still lacking is a physical picture of how the large, wrapping loops of current disorder 
the vacuum and produce the string tension. The fact that these loops extend over the 
whole lattice suggest that there are low mass (perhaps massless) magnetically charged 
excitations present in the confined phase. We plan to report elsewhere on this question 
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as well as how the magnetic current screens itself. The results obtained in our work are 
likely to have an impact on the monopole approach to confinement in non-Abelian gauge 
theories. 
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Figure 1: The quark potential calculated using the magnetic current configurations 
for = 1.0103 (triangles) and jS^y = 1.0114 (squares). The solid lines are lin- 
ear-plus- Coulomb fits to the potentials. 
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Figure 2: The potential calculated using only the wrapping monopole loops (triangles) 
and using only the simple monopole loops (squares). The photon contribution from Wph 
has not been included. 



13 



Figure 3: Comparison of the quark potential calculated using the full magnetic current 
and photons (triangles) with the potential obtained by summing the potentials determined 
separately from photons, wrapped loops, and simple loops (squares). The solid line is the 
hnear-plus-Coulomb fit from Fig. (1) for Pw — 1.0103. 
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